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1.  INTRODUCTION 

When  multilayer  interference  filters  are  used  in  a conical  light  beam  the 
effective  spectral  reflectance  and  transmission  characteristics  differ  from 
those  which  apply  when  only  parallel  pencils  at  normal  incidence  are  involved. 

Some  years  ago  a computer  programfref . 1)  was  written  to  analyse  the  effect 
of  conical  light  beams  on  the  performance  of  narrow  passband  filters.  In  this 
program  it  was  assumed  that  the  filter  was  illuminated  through  a variable 
circular  aperture  for  various  conditions  of  filter  orientation  within  the  beam, 
so  that  the  effect  of  different  cone  angles  and  filter  tilt  angles  could  be 
studied.  This  situation  occurs  when  a filter  is  placed  in  the  converging  beam 
from  a lens  system  which  is  focusing  a point  source  object.  The  results  obtained 
from  this  program  can  be  used  to  determine  the  effect  of  cone  angle  and  tilt  on 
the  location  and  spectral  width  of  the  filter  passband. 

When  a multilayer  filter  is  incorporated  into  an  optical  system  with  an 
extended  field  of  view  different  sections  of  the  object  plane  give  rise  to  conical 
beams  that  are  tilted  at  different  angles  to  the  surface  of  the  filter.  This  is 
a more  complicated  situation  than  was  considered  previously  and  the  filter  per- 
formance will  depend  on  the  integrated  effect  of  all  the  tilted  conical  beams 
incident  upon  it. 

It  recently  became  necessary  to  design  a narrow  passband  filter  for  use  in  a 
radiometer  system  having  a defined  aperture  and  field  of  view.  The  spectral 
width  and  location  of  the  filter  passband  when  used  in  the  lens  system  was  of 
some  importance  and  was  required  to  lie  within  specified  limits.  In  order  to 
establish  the  design  restrictions  for  the  optical  system,  some  additions  and 
alterations  to  the  computer  program  of  reference  1 were  made  to  include  the 
effect  of  a defined  field  of  view. 

This  report  describes  the  modifications  made  to  the  previous  computer  program 
and  discusses  the  effect  of  variations  of  aperture  and  field  of  view  on  the  per- 
formance of  all-dieletric  Fabry-Perot  interference  filters  when  they  are  used 
in  radiometers  and  other  similar  optical  systems. 


2.  DESCRIPTION  OF  THE  COMPUTATION 


2 . 1 General 

A brief  discussion  on  the  optical  arrangement  used  in  simple  filter 
radiometers  will  illustrate  the  analysis  problem.  It  is  assumed  in  this 
analysis  that  the  spectral  radiance  of  the  field  of  view  is  uniform  or 
randomly  distributed  such  that  the  average  spectral  radiance  is  uniform 
on  a course  macroscopic  scale.  The  simplest  use  of  a lens,  filter, 
detector  combination  is  shown  in  figure  1.  The  lens  produces  a real  image 
of  the  field  of  view  at  the  detector  surface,  and  the  filter  is  assumed  to 
lie  in  any  plane  normal  to  the  axis  between  the  lens  and  the  detector. 

The  image  formed  is  a faithful  reproduction  of  a distant  ooject  and  will 
contain  the  small  intensity  variations  of  the  object  field.  This  is  un- 
desirable in  a radiometer  because  there  may  be  variations  in  sensitivity 
across  the  detector.  It  is  preferable  to  irradiate  the  detector  surface 
uniformly  and  this  is  done  in  figure  2 by  the  use  of  a field  lens  which 
images  the  objective  aperture  onto  the  detector.  The  light  entering  the 
objective  lens  comes  from  all  parts  of  the  field  of  view  and  it  is  therefore 
uniformly  illuminated  across  the  aperture.  The  objective  lens  of  aperture 
"A"  produces  an  image  of  the  object  at  the  position  occupied  by  the  field 
lens.  The  focal  length  and  position  of  the  field  lens  is  such  that  it 
produces  an  appropriately  sized  image  of  the  aperture  "A"  at  the  detector. 
Although  this  optical  system  contains  two  lenses,  for  the  purpose  of 
analysing  filter  performance  it  can  be  considered  as  equivalent  to  a single 
lens  with  an  aperture  of  diameter  (P  < D)  located  at  a distance  Q from  the 
detector  of  diameter  d.  The  important  parameters  which  affect  filter 
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performance  are  the  half  cone  angle  a = tan'1  — ■ 

d/2  ^ 

has  a maximum  value  of  t = tan-1  (-^-). 


and  the  tilt  angle  which 


The  geometry  of  the  equivalent  optical  system  to  be  considered  is  that 
shown  in  the  lower  drawing  of  figure  2.  At  the  centre  of  the  detector  the 
irradiation  is  by  a convergent  cone  of  light  of  half-angle  a given  by 
tan-1  (P/2),  with  the  axis  of  the  cone  at  normal  incidence  to  the  detector 


Q 


surface.  At  the  edge  of  the  detector  the  irradiation  is  by  a cone  whose 
axis  is  tilted  at  an  angle  t given  by  tan  1 (d/2) . At  other  points  on  the 

Q 

detector  surface  there  is  some  combination  of  the  cone  angle  and  a tilt 
angle,  the  former  decreasing  and  the  latter  increasing  with  the  distance 
from  the  detector  centre. 

The  object  of  the  computer  program  is  to  calculate  the  radiation  falling 
on  the  detector  (and  also  the  filter)  as  a function  of  angle  of  incidence 
and  hence  to  calculate  the  effective  spectral  transmission  of  the  filter. 

This  is  done  in  two  steps.  Firstly,  the  detector  surface  is  divided  into 
annular  rings  centred  on  the  axis  of  the  optical  system  and  having  equal 
increments  of  radius.  For  each  ring  a mean  cone  and  tilt  angle  is  computed 
for  the  geometry  of  the  lens  system  under  consideration.  Having  established 
the  mean  cone  and  tilt  angle  associated  with  each  ring  the  original  program 
(ref.l)  is  used  to  determine  the  distribution  of  the  radiation  over  the 
incidence  angles  contained  within  each  tilted  cone.  This  population  is 
given  a weighting  factor  proportional  to  the  area  of  the  ring  on  the  detector 
surface  with  which  the  cone  is  associated.  The  analysis  is  performed  for 
each  ring  in  turn  from  the  centre  of  the  detector  to  the  edge.  The  angular 
distribution  of  radiation  for  all  rings  is  summed  so  that  the  final  output 
from  the  first  part  of  the  program  is  an  array  of  angles  and  the  relative 
amount  of  radiation  incident  on  the  detector,  and  also  the  filter,  at  those 
angles. 

The  second  half  of  the  program  computes  the  spectral  transmission  of  the 
interference  filter  at  each  of  the  incidence  angles  specified  by  the  first 
half.  These  spectral  transmission  curves  are  thru  weighted  by  the  angular 
distribution  of  the  incident  radiation  and  summed  to  give  the  effective 
spectral  transmission  curve  of  the  filter  when  operated  in  the  lens  system. 
The  steps  in  the  computation  will  now  be  discussed  in  greater  detail. 


2.2  Determination  of  the  angles  of  incidence 


The  geometry  of  the  lens  system  is  defined  for  the  program  by  the  equiv- 
alent aperture  P,  the  separation  Q and  the  detector  diameter  d.  Referring 
to  the  computer  program  flow  chart  shown  in  figure  3 the  first  step  in  the 
calculation  is  to  determine  the  angles  at  which  the  spectral  performance 
of  the  filter  will  be  calculated  later  in  the  program.  This  is  done  by 
computing  the  maximum  angle  of  incidence  of  the  light  falling  on  the  filter 
surface  and  dividing  this  into  ten  incremental  angles.  Referring  to 
figure  2 the  maximum  angle  of  incidence  of  the  cone  is: 


0 max  = tan'1 


rd/2  + P/2. 

1 Q } 


The  minimum  angle  is  normally  zero  for  concentric  lens  systems.  Ten  angles 
have  been  found  to  give  satisfactory  results  from  the  calculations.  Using 
smaller  increments  of  angle  does  not  change  the  result  significantly. 
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The  next  program  step  divides  the  detector  surface  into  concentric  rings 
and  calculates  the  cone  and  tilt  angles  which  are  associated  with  each  ring. 
Let  the  mean  radius  of  the  ring  be  r.  Then  from  figure  4 the  tilt  angle 
t = tan-1  (r/Q). 

The  cone  angle  is  obtained  by  calculating  the  maximum  and  minimum  angles 
in  the  cone.  Then  the  half  cone  angle,  which  is  half  the  included  angle 
(8  - i?)  is: 


cone  angle  =0.5  90  + tan 


rP/2 


- tan 


r + 


p/j) J when 


< P/2 


cone  angle  =0.5 


tan 


tan'1  (frS/2') ) when  r > P/2 


The  remaining  calculations  were  fully  described  in  reference  1,  but 
because  of  changes  introduced  in  the  notation  an  account  of  the  steps 
involved  will  be  given  here. 

The  angle  of  incidence  at  the  filter  is  then  evaluated  for  rays  from  a 
number  of  regularly  spaced  locations  within  the  aperture  to  establish  the 
angular  distribution  of  incident  radiation  within  the  range  of  incidence 
angles. 

Figure  5 shows  the  geometry  of  the  system  under  consideration.  The 
upper  diagram  shows  a conical  beam  of  light  emanating  from  a circular 
aperture  of  radius  D/2  in  a plane  perpendicular  to  the  cone  axis  and 
focusing  at  the  point  0 on  the  filter.  One  half  of  the  plane  is  shown 
divided  into  an  array  of  squares  located  at  positions  specified  by  the 
coordinate  system,  x,y.  When  considering  the  case  for  that  portion  of 
the  filter  at  normal  incidence  to  the  axis  of  the  irradiating  cone,  the 
aperture  is  assumed  to  be  in  the  plane  as  shown  in  the  top  diagram.  To 
consider  those  parts  of  the  filter  tilted  to  the  axis  of  the  irradiating 
cone  it  is  assumed  that  the  aperture  ar.d  filter  are  rotated  about  the 
x-axis  by  the  angle  of  tilt.  The  top  half  of  the  aperture  will  then  lie 
behind  the  plane  containing  the  grid  reference  x,y  and  the  lower  half  will 
lie  in  front  of  it.  To  establish  the  distribution  of  incidence  angles 
for  rays  passing  through  all  parts  of  the  aperture  the  incidence  angle  is 
calculated  at  each  position  on  the  aperture  where  it  is  cut  by  lines  from 
the  origin  through  the  grid  positions  x,y. 

All  the  squares  in  the  array  which  lie  within  the  circular  aperture  are 
included  in  the  analysis  i.e.  at  normal  incidence  all  positions  where 


x2  + y2  < P/2 


For  tilted  cones  (the  general  case)  the  grid  positions  are  projected 
back  to  the  vertical  to  test  if  they  lie  within  the  circle  of  radius  P/2. 
They  are  included  when: 


and 


yx2  + (y.cos  t)2  < P/2  for  +VE  y 


Vx2  + (y/cos  t)2  < P/2  for  -VE  y 
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The  lower  diagram  on  figure  5 shows  an  enlarged  view  of  ABCD  a portion 
of  the  surface  of  the  grid  reference  plane  with  the  aperture  tilted  about 
the  x-axis,  which  is  the  most  general  case. 

The  incidence  angle  9 which  is  to  be  calculated  is  the  angle  between  the 
incident  light  and  the  normal  at  the  filter  surface.  For  light  passing 
through  the  aperture  at  the  position  C',  the  angle  0 is  calculated  as 
follows: 

From  figure  5 the  vertical  angle: 

P = tan-1  (y/OA) 

The  angle  is: 


^ - tan'1  (x/(OA  sec  j3)) 


From  the  triangle  OC’F’ 


0 = cos  1 OF* 
OC' 


Where  the  normal  to  the  filter  and  aperture  planes 


OF'  = OB'  cos  (P  + t)  from  triangle  OB'F' 


and 


OC'  = OB'. sec  from  triangle  OB'C' 


The  incidence  angle  on  the  filter  is  then: 

9 = cos-1  (cosQ?  + t)cos  (i/0) 

A similar  procedure  is  used  to  calculate  the  incidence  angle  from  the  position 
(x,-y)  below  the  axis. 

As  the  incidence  angle  associated  with  each  grid  position  is  evaluated, 
the  nearest  assigned  angle  is  found  and  an  integer  is  added  to  a store 
associated  with  it.  When  the  aperture  area  has  been  completely  scanned, 
the  population  held  in  the  store  associated  with  each  assigned  incidence 
angle  gives  the  number  of  grid  positions  where  this  assigned  value  of 
incidence  angle  was  encountered.  For  the  calculations  discussed  in  this 
paper,  x and  y were  arbitrarily  given  maximum  values  of  ±50  (i.e.  2500  ir 
unit  squares  in  the  irradiated  area) . 

This  analysis  is  carried  out  for  cone  and  tilt  angles  specified  for 
each  of  the  rings  on  the  filter  surface.  The  total  population  at  each 
assigned  angle  after  the  completion  of  this  analysis  is  multiplied  by  a 
factor  proportional  to  the  area  of  each  ring,  and  normalised  to  give  the 
angular  distribution  of  the  irradiance  incident  on  the  filter. 

Having  established  a range  of  incidence  angles  and  their  associated 
relative  intensity  figures  for  light  incident  on  the  filter,  the  next  step 
is  to  calculate  the  performance  of  the  filter  at  each  of  these  incidence 
angles,  multiply  the  values  of  the  reflection  or  transmission  coefficient 
by  the  angular  distribution,  and  sum  the  ten  performance  curves  to  give  the 
integrated  effect  for  the  irradiation  conditions  over  the  wavelength  range 
of  interest. 
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The  methods  used  to  calculate  the  reflection  and  transmission  coeffi- 
cients of  the  filter  were  described  fully  in  reference  2,  but  for 
completeness  a brief  description  is  given  in  the  following  section. 

2.3  Computation  of  the  filter  characteristics 

The  method  of  calculation  and  the  notation  used  are  those  suggested 
by  Weinstein (ref . 3)  and  are  illustrated  in  figure  6. 

The  tangential  components  of  the  electric  and  magnetic  field  strength 
vectors  E(J+1)  and  H(J+1)  immediately  to  the  left  of  the  Jth  layer  of  the 
thin  film  system  may  be  expressed  in  matrix  from  as: 


/ E(J+1)  \ / cos  g(J)  _i sin  g(J)  \ ? E(J)\ 

\ H(J+1)  / \ i u(J)  sin  g(J)  U(J)  cos  g(J)  / ' \ H(J)/ 

or 

( E<J*»  \ = H(J)  . ( 

\ H(J.l)  j V H(J)/ 


Where  M(J)  is  the  characteristic  matrix  of  the  Jth  layer. 

u(J)  may  be  considered  as  an  effective  refractive  index  which 

varies  with  angle  of  incidence  and  the  state  of  polarisation 

g(J)  is  the  effective  optical  thickness  of  the  layer,  and  is  given  by: 
g(J)  = 2»/x  . n(J)  h(J)  cos  fl(J) 

X is  the  wavelength  of  the  radiation  under  consideration 

(nanometres) 

n(J)  is  the  refractive  index  of  the  layer  material 

h(J)  is  the  thickness  of  the  layer  (nanometres) 

0(J)  is  the  angle  of  incidence  of  the  radiation  in  the  Jth  layer 

By  successively  applying  this  relationship  to  each  layer  of  a thin  film 
construction,  an  expression  is  obtained  for  the  values  of  E(P)  and  H(P), 
the  field  strength  immediately  to  the  left  of  the  outside  layer  of  the 
filter,  (figure  6)  in  terms  of  the  matrices  of  the  layers  and  E(l)  and 
H(l). 

Thus 


E(p)  \ ! E(l)\ 

HP)  ) - M^-15  ' • V H(1V 


or  in  general 


E(P) 

H(P) 


= M(product)  . H(1^ 
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Where  M(product)  is  the  matrix  formed  by  multiplying  the  characteristic 
matrix  of  each  layer  by  the  next  succeeding  layer  matrix  and  continuing  the 
process  through  the  filter  from  layer  1 to  layer  (p-1),  the  final  layer. 

Having  calculated  the  product  matrix  of  the  complete  thin  film  system, 
the  values  of  E(p)  and  H(p)  may  be  calculated  by  using  the  relationship 


H 


(J) 


(J) 


'(J) 


and  assuming  the  initial  conditions 


ECl)  - S(s)  ■ 1 


H(l)  = n(s)  • E(s)  * "(,) 


and  evaluating  the  matrix.  (The  subscript  s denotes  the  value  of  E and  n 
pertaining  to  the  substrate.) 

If  absorbing  materials  are  used  in  the  filter,  the  complex  form  of  the 
refractive  index  n = N -iK  must  be  used,  where  K is  the  absorption 

4T  V 

coefficient,  which  leads  to  a reduction  of  the  complex  amplitude  by  e 2 
per  vacuum  wavelength  traversed  in  the  direction  of  propagation.  The 
complex  refractive  index  may  be  substituted  directly  into  the  above 
equations  instead  of  n(J),  and  the  values  of  E(p)  and  H(p)  which  are 
produced  will  then  be  complex.  See  for  instance  reference  4. 

The  reflection  and  transmission  coefficients  may  then  be  calculated 
(ref. 5)  from: 


_ E(p)  - H(p)  / N(p)  2 
E (p)  + H(p)  / N(p) 


t*  — - V.S  J 2 

' N(p)  )E(p)  + H(p)/N(p)| 


A flow  chart  for  the  matrix  calculation  and  graph  plotting  is  shown  in 
figure  7. 

Because  of  the  large  number  of  computations  performs*  by  the  program  in 
any  single  run,  output  from  it  is  normally  restricted  to  the  final  graph 
showing  the  effective  spectral  transmission  of  the  filter  under  the  specified 
illumination  conditions.  The  optical  constants  of  the  materials  used  in 
the  thin  film  designs  under  analysis  are  read  from  a data  set  (described 
in  reference  6)  which  contains  a library  of  the  optical  constants  of  the 
substrate  and  evaporation  materials  used  in  the  manufacture  of  thin  film 
interference  filters.  The  library  is  updated  when  new  materials  are 
used  or  the  evaporation  conditions  are  changed,  as  these  effect  the  optical 
constants  of  materials  in  thin  film  form.  The  filter  description,  wave- 
length and  plotting  specifications  and  the  optical  system  geometry  are  read 
from  a separate  data  set  allocated  by  the  controlling  program. 
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As  an  illustration  of  the  type  of  result  obtained  from  the  computation 
two  examples  of  narrow  passband  filter  performance  in  focused  lens 
systems  are  given  in  the  next  section. 


3.  EXAMPLES  OF  CHANGES  IN  FILTER  PERFORMANCE 

3.1  Description  of  the  filter 

The  filter  construction  under  consideration  is  an  all -dielectric 
Fabry-Perot  design  which  has  a layer  system  as  follows: 


substrate/HLHLHLHL. 12H. LHLHLHLH/air 


where  H and  L represent  high  and  low  refractive  index  layers  with  an  optical 
thickness  of  A/4  at  1060  nm.  At  normal  incidence  the  filter  has  a half 
width  (T  = 50%)  of  1.3  nm  and  when  the  substrate  back  surface  reflection 
is  omitted  it  has  a computed  transmission  of  96.6%  at  1060  nm.  The  normal 
incidence  spectral  transmission  curve  for  the  filter  passband  is  shown  as 
the  highest  curve  in  figure  10. 

The  two  cases  considered  in  this  section  are  variation  of  performance 
due  to  changes  in  aperture  (i.e.  variations  in  the  F-number)  and  due  to 
changes  in  field  angle  (diameter  of  detector).  The  filter  can  be 
considered  to  lie  in  any  plane  normal  to  the  axis  between  the  aperture 
and  the  detector. 

3.2  Effect  of  variation  of  aperture 

Figure  8 shows  a series  of  spectral  transmission  curves  for  the  filter 
when  the  aperture  of  an  irradiation  system  with  a constant  field  angle  is 
changed.  The  dimensions  chosen  are  a variable  aperture  of  20  mm  maximum 
diameter,  a fixed  detector  diameter  of  10  mm  and  a separation  of  80  mm, 
giving  an  F/4  system  at  full  aperture  and  a semi-field  angle  of  3.58°  at 
full  aperture.  Under  these  conditions,  the  filter  peak  transmission  drops 
to  52%  and  shifts  to  1059  nm,  and  the  half  width  increases  to  3.0  nm. 

As  the  aperture  is  decreased  the  filter  performance  improves.  At  an 
aperture  of  10  mm  CP/8)  the  transmission  rises  to  76%  at  1059.6  nm  and  the 
half  width  is  1.9  nm.  At  an  aperture  of  5 mm  (F/16)  the  peak  transmission 
is  87%  at  1059.7  nm  and  the  half  width  is  1.7  nm. 

3.3  Effect  of  change  in  field  angle 

The  effect  of  field  angle  changes  was  calculated  for  a 10  mm  diameter 
aperture,  a variable  filter  diameter  and  a separation  of  80  mm  i.e.  an 
F/8  system.  These  parameters  were  chosen  to  give  a reasonable  filter 
performance  as  shown  by  figure  8.  Figure  9 shows  that  in  the  worst  case 

considered,  a 20  mm  diameter  detector  (semi-field  angle  of  7.12°)  the  peak 
transmission  is  55%  at  1059.2  nm  and  the  half  width  is  2.9  nm.  The 
smallest  detector  diameter,  5 mm  (semi-field  angle  1.79°)  gives  a peak 
transmission  of  85%  at  1059.7  nm  and  a half  width  of  1.5  nm. 

Figure  10  shows  the  effect  of  varying  the  F-number  of  the  optical  system 
when  the  aperture  and  detector  are  the  same  diameter.  The  trends  are 
similar  to  those  shown  in  figure  7 for  variations  in  aperture  with  a fixed 
detector  diameter.  Providing  the  filter  is  used  at  an  aperture  no  wider 
than  F/8  the  performance  is  not  markedly  degraded. 
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3.4  General  comments 

The  parameters  which  were  varied  in  Section  3.2  and  3.3  have  the  same 
general  effect  on  the  shape  and  location  of  the  passband.  As  the  aperture 
of  the  system  is  increased  the  bandwidth  increases,  the  transmission  drops 
and  the  passband  shifts  to  a shorter  wavelength.  The  reason  for  these 
changes  is  simply  that  the  filter  is  being  irradiated  with  light  at  a 
variety  of  different  incident  angles  and,  as  predicted  by  theory  the  filter 
passband  shifts  to  a shorter  wavelength  as  the  angle  of  incidence  of  the 
light  increases. 

Several  of  the  transmission  curves  shown  on  figures  8 and  9 represent 
cases  where  the  aperture  size  and  detector  size  have  been  interchanged. 

From  consideration  of  the  relevant  ray  diagrams  one  would  expect  a similar 
filter  performance  in  these  cases.  The  small  differences  that  are 
apparent  between  the  curves  for  interchanged  dimensions  are  due  to  rounding 
errors  in  the  numerical  integration  process.  Taking  smaller  increments 
of  angles  and  radii  on  the  detector  surface  would  reduce  these  errors  at 
the  expense  of  greater  computer  storage  requirements  during  the  calculation. 
These  errors  are  smaller  than  the  differences  which  often  occur  between 
computed  and  measured  filter  performance  due  to  inaccuracies  of  manufacture 
and  errors  in  the  refractive  indices  used  in  the  computations. 

For  the  particular  case  of  figure  8 the  angular  energy  distribution  of 
the  incident  radiation  has  been  plotted  for  each  of  the  transmission  curves 
shown.  Figure  11  shows  the  normalised  angular  energy  distribution  which 
is  the  output  from  the  first  half  of  the  computation.  The  characteristic 
shape  of  the  distribution  is  the  same  in  each  case,  but  the  angular  spread 
increases  with  the  aperture  of  the  lens  system.  For  the  smallest  aperture 
the  transmission  peaks  calculated  at  each  of  the  incidence  angles  shown  in 
curve  A (figure  11)  overlap  to  produce  a narrow,  high  transmission  region. 

As  the  aperture  is  increased  the  range  of  incidence  angles  is  greater  and 
the  effective  transmission  peaks  are  spread  over  a greater  wavelength  range. 
The  resultant  passband  is  then  wider  and  lower  in  transmission  because  there 
is  less  overlap  in  the  individual  transmission  peaks. 

The  energy  distribution  over  the  range  of  incidence  angles  changes  with 
the  geometry  of  the  optical  system.  This  is  demonstrated  in  figure  12 
which  shows  the  output  from  the  first  part  of  the  computation  for  a variety 
of  different  ratios  of  aperture  to  detector  diameter  at  a fixed  separation. 

The  curves  have  been  normalised  to  remove  the  effect  of  different  detector 
diameters  and  they  are  plotted  against  a normalised  range  of  incidence  angles. 
When  the  aperture  and  detector  diameters  are  equal  the  distribution  is 
approximately  symmetrical  over  the  range  of  angles  as  shown  by  curve  A. 

When  the  aperture  and  detector  diameters  are  different  there  is  more  energy 
incident  at  higher  angles  of  incidence  as  shown  by  curves  B and  C. 

If  a filter  is  to  be  used  under  these  irradiation  conditions  it  would  be 
necessary  to  make  it  at  a longer  wavelength  them  required  so  that  the  wave- 
length shift  of  the  transmission  peak  due  to  the  convergence  of  the  lens 
system  moves  it  to  the  correct  wavelength.  A compromise  would  be  necessary 
if  a variable  aperture  is  involved  because  the  location  of  the  peak  shifts 
as  the  aperture  is  changed. 


4.  CONCLUSION 

This  computation  has  been  applied  to  filters  intended  for  use  in  a variety  of 
optical  systems  and  has  provided  a useful  method  of  assessing  the  effective 
filter  passband  location  and  bandwidth.  The  results  of  the  computation  have 
been  useful  in  specifying  the  optimum  control  wavelength  for  the  manufacture  of 
narrow  passband  filters  for  use  in  optical  systems  having  an  extended  field  of 
view  over  which  the  spectral  radiance  can  be  considered  uniform  or  to  be  randomly 
distributed  such  that  the  average  spectral  radiance  is  uniform  on  a course 
macroscopic  scale. 
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NOTATION 

field  lens  diameter 

tangential  component  of  electric  field  strength 
tangential  component  of  magnetic  field  strength 
incidence  angle 

integer  denoting  layer  numbers 
absorption  coefficient 
a two  by  two  matrix 

real  part  of  the  complex  refractive  index 

diameter  of  aperture 

separation  of  aperture  and  detector 

reflection  coefficient 

transmission  coefficient 

filter  diameter 

focal  length 

effective  optical  thickness  of  a layer 

thickness  of  a layer 

complex  constant 

refractive  index  of  a layer 

mean  radius  of  ring 

subscript  denoting  the  substrate 

angle  of  tilt 

effective  refractive  index 

array  co-ordinates 

half  cone  angle 
vertical  angle 
vertical  angle 
vertical  angle 


V 
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“ angle  of  incidence 

^ wavelength 

P layer  number 

^ horizontal  angle 
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Figure  1.  Simple  radiometer  lens  system 


Figure  2.  Two  lens  radiometer  and  equivalent  diagram 
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Figure  7.  Flow  chart  of  the  filter  transmission  calculation 
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Figure  10.  Filter  performance  changes  due  to  variation  in  F number 
of  optical  system 
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